RADIO SCIENCE Journal of Research NBS/USNC-URSI 
Vol. 69D, No. 6, June 1965 

Propagation of Waves Across a Magnetoplasma- 

Vacuum Boundary ' 

Robert L. Gailawa 
Central Radio Propagation Laboratory, National Bureau of Standards, Boulder, Colo. 

(Received January 13, 1965; revised January 23, 1965) 

The coupling of electroacoustic and electromagnetic energy at a magnetoplasma- vacuum boundary 
is studied via the plane wave solution of Maxwell's equations. The effect of the static magnetic field 
is illustrated graphically. It is found that the coupling depends significantly on the angles, increased 
coupling being in favor of the grazing angles of transmission. The electron density of the plasma is 
seen to be an important parameter in most cases. 



1. Introduction 

In recent years there has been a considerable amount of energy expended in an effort to 
understand more completely the coupling mechanism which governs the conversion of electro- 
magnetic energy at a plasma density discontinuity. In this general problem is included the 
interaction of electromagnetic energy with the plasma and a completely general account of the 
mechanism is very difficult. In order to reduce the problem to a tractable form, simplifying 
(nonphysical) assumptions are made which necessarily restrict the domain of application. Thus, 
for example, the linearized form of the hydrodynamic equations are used [Oster, I960]; this restricts 
the amplitudes to small (perturbation) values. 

The problem of energy conversion due to a boundary has received some attention by con- 
sidering the reflection from and transmission into an isotropic plasma half-space [Hessel et al, 
1962; Kritz and Mintzer, 1960; and Wait, 1964a], and more recently by Wait [1964b] for a magneto- 
plasma half-space. In these cases the problem has been one of characterizing the wave in a plasma 
half-space due to a wave incident from a vacuum. In recent years there has been an effort to 
understand the reciprocal problem, that of determining the consequences of a wave incident from 
the plasma onto a plasma-vacuum boundary [Field, 1956; Johler, 1964]. 

2. Formulation of the Equations 

In this study a plane interface is postulated which separates a half-space of plasma and a half- 
space of vacuum. The plasma is taken to be a one-component, uniform, electron fluid; i.e., heavy 
ionic motion is neglected. However, collisions between electrons and neutral particles are 
accounted for by a constant collision frequency, v. In addition, the finite compressibility is 
considered. The plasma is anisotropic by virtue of a steady magnetic field which permeates the 
region. Effort will be directed towards a description of the conversion factor which characterizes 
the (transverse) electromagnetic field existing in the vacuum due to a longitudinal (electroacoustic) 
wave from the plasma onto the interface. The wave magnitudes are assumed to be very small, 
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permitting the use of linearized hydrodynamic equations [Oster, 1960] which are given by 

mre (v + |-jF=noe(E+FxF)-Vp, (1) 

u*mnoV-V=-&, (2) 

Vx£=- Mo f, (3) 

VX// = €o ^ + no6^ (4) 

where e and tti are the charge and mass of the electron, respectively; /xo and 6 are the permeability 

and dielectric constant of free space, respectively; u is the speed of sound in the electron gas, 

KT 

u 2 = , X = Boltzmann's constant, T= absolute temperature; p is the deviation of the electron 

m 

pressure from the mean and V is the mean electron velocity; /i is the constant electron number 
density. E and H are, as usual, the electric and magnetic fields. B is the steady magnetic 
induction. Without loss of generality, a (surpressed) time factor exp (icot) is assumed. 



3. Homogeneous, Infinite Plasma 

In order to understand the propagation mechanism, it is appropriate to seek a plane-wave 
solution which is proportional to exp (— ikn • r), where n is the unit vector in the direction of propa- 
gation and r is the vector to the variable point. Under these conditions, 



V =—ikn. 



(5) 



Unit vectors a and 1 are also introduced such that n, a, andTform a right-handed triad and the basis 
of a rectangular coordinate system. Attention will be restricted to the case where B = aB . 
Then for the assumed time variation 



c 2 + F^ 



+ [nEn + aE a + tE t ] [k 2 - k 2 ] 



+ tE n 



or 



where 



,[- 

iucjc , co c k 2 1 - [ iiococ , Q) c k 2 1 ~ r (*>ck 2 

— — I — —T-\-nEt — — H — — — \-tEn— 7-r-=0, (6) 

E a [k 2 -K 2 ] = 0, (7b) 

£ 4 z ^ £ ]"^ [ ^ 2 ~ ,c2]=o ' (7c) 



2 2 . CO/Xo€oCog 

k 2 = a>Vo€ - i , : p , (8) 

2 V+JO) 

v 2 = ~ , (9) 

l(OjJLo€o 

808 



10 



2 =- 
P 



a>c=—eBolm, (11) 

and 

E = nEn + aE a +JE t . (12) 

It is clear from (7) that E a is not coupled to the remaining components and has a propagation 
constant 

k 2 = K 2 . 

The field components E n and E t are coupled, however, by virtue of the static magnetic field. The 
modes which propagate are then quasi-longitudinal and quasi-transverse. This description 
indicates the fact that neither a purely longitudinal wave nor a purely transverse wave can be 
supported in this medium. If the static magnetic field vanishes, it is readily seen that the medium 
can support three waves, each propagating independently. The wave which is polarized in the 
direction of propagation (the longitudinal wave) has propagation constant 

v 2 

k 2 = K 2 ~ 
UT 

The transverse waves, under these conditions, each has propagation constant 

k 2 =K 2 . 

It is interesting to note that there is no thermal correction for the wave number, A:, for the 
mode which is polarized in the direction of the magnetostatic field. However, there is a thermal 
correction for the wave numbers associated with the other two waves by virtue of the term u 2 \v 2 
which appears in (7a). If the magnetostatic field vanishes, the thermal correction for the com- 
ponent E t vanishes and the index of refraction for the transverse waves is the same as it would be 
for a cold (incompressible) plasma. 



4. Case of Zero Static Magnetic Field 

In order to develop concepts which will be useful in later work, consideration is given to the 
case of vanishing static magnetic field. In that case, (1) is modified to the extent that B is set 
equal to zero. Following earlier workers [Hessel and Shmoys, 1961; Hessel, Marcuvitz, and 
Shmoys, 1962; Chen, 1964; Field, 1956; Yildiz and Silberg, 1964], it is then convenient to split 
E into longitudinal and transverse components. Let 

E = El H~ Et, 
where 

(13a) 
(13b) 



Vx£ t = 


= 0, 


V 


■E r = 
809 


= 0. 



Then by (2), (4), and (13b), 



V-E L = 



L—~Z p. 

U*€oTTl 



In addition, 



V xE T = —i(DfJL()H. 



(14a) 



(14b) 



From these equations, it is evident that El is in the direction of propagation (and hence is longitu- 
dinal) while Ej is entirely transverse to the direction of propagation. Equations (13) and (14) com- 
pletely specify E L and Et. The equation satisfied by E is found in the usual manner by taking the 
curl of (3): 

- V 2 E T -k 2 ^T-k^ L + io)fi noeV= 0, (15) 

where k\ = co 2 /x €o. The velocity Fmay also be broken into longitudinal and transverse components: 



V=V L +V T . 
Then from (1) 

mno(v -\- ia))VL + mno(v + ia))VT = ti^El -f rioeEr — Vp. 
It is convenient and congruous to let 

mno(v-\- ii^)Vh — n^eEi — Vp, 
mn {v 4- iw)Vr — n^eEr. 
Then (15) becomes 



[-V 2 -k 2 ]Et~- 
Taking the curl of (18) yields 



k 2 -^ 



E L + : \_. \ ^P- 



[-V 2 "K 2 ]VX£ T =0. 
Since the curl of Et is nonzero unless Et is zero, this demands that 



V 2 £r + ^ 



\l !-L-l 



£V = 0. 



Taking the divergence of (18) yields 



ur 



CO Z CO 



p = 0. 



(16) 

(17a) 
(17b) 

(18) 



(19) 



(20) 



Equations (19) and (20) agree with the equations of Wait [1964a]. They show explicitly the propa- 
gation constants associated with the transverse electromagnetic wave and the pressure wave. 
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Furthermore, the pressure wave has an associated longitudinal electromagnetic wave by virtue 
of (14a). The two types of waves travel independently and with propagation constants 



K L~ 2 





l-l 


o> 2 


. V 

I — 



cot 



co(v + ito)_ 

n 2 



U Z CO u l 



(21a) 
(21b) 



It is interesting to note that in the absence of collisions, the ratio of the two propagation constants 
is just the ratio c/u. 



5. Finite Static Magnetic Field 

The set of equations (7a) and (7b) leads to solutions for k by equating the determinant of the 
coefficients to zero. This leads to 



k2+k2 ^ — ^L 



F = - 



v 2 v x (v + i(o) 






V\ p« l v\v + io>)) % 2 \ v* 



(22) 



v- 



It can readily be shown that in the limit of zero static magnetic field, the upper sign leads to 

k2 = l < 2 ^ (23a) 



and the lower sign yields 



k 2 = K 2 ' 



(23b) 



Comparison of (23) with (21) indicates that the upper sign should be associated with the quasi- 
longitudinal wave (ki) and the lower sign with the quasi-transverse wave (At). 



6. Coupling at a Plasma- Vacuum Boundary 

The mechanism of coupling at a boundary can be studied by examining the mechanism of 
reflection and transmission at a plasma-vacuum boundary. Attention is restricted to the case of 
coupling from an electroacoustic wave. In order to do this, a quasi-longitudinal plane wave is 
assumed to be incident on the boundary from the plasma (see fig. 1). Such a wave could originate 
in the solar corona as discussed by Field [1956]. 

The magnetostatic field is taken to be in the y-direction, as discussed. 2 It is supposed that 
the incident wave is given by 

ELi = Ew\a x cos 6 + ciz sin d~\~R{—a x sin + a z cos 6)] • exp [—ikiix cos + z sin 0)], (24) 



2 The treatment here is abbreviated since a closely related problem has already been considered by Wait [1964b]. 
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MEDIUM 1 
£ , u o , n o 




MEDIUM 2 



Figure 1. Coordinate system showing the two media and 
illustrating the angles of incidence, reflection, and 
transmission. 



/here 



The reflected longitudinal wave is 



R 



lOJCJc 



V 2 (kf-K 2 ) 



(25) 



ELr = Eir\—a x cos 6 + a z sin + R{—a x sin 6 — a z cos 0)] • exp [—ik L (—x cos d + z sin 0)]. 

The reflected transverse wave can be of two types: 

Ety = E' Tr \— a x sin (p\ — a z cos <p\ + T(— a x cos (pi + a* sin (pi)] 

• exp [— ikn(—x cos <p\-\~z sin <pi)] . 
E { p r = a y E" Tr exp [— ikn(— x cos <pi+z sin <£i)], 
where 

— icococ 
The transmitted wave can likewise be of two types: 

Efy = E r TT [a x sin (p2~a z cos ^2] exp [— ikrz{x cos (^2 + 2 sin <£>2)], 
E^ = a y Ej T exp [— ikrzix cos <^?2 H~ -2: sin (^2)]. 



(26) 



(27) 
(28) 

(29) 

(30) 
(31) 



The notation used on the angles is as follows: 6 is used in association with the longitudinal waves and 
<p is used in specifying the transverse waves (see fig. 1). The subscript 1 or 2 is used to designate 
plasma or vacuum, respectively. 
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In order to determine the magnitudes of the unknown amplitudes, the boundary conditions are 
applied at the interface. Because the interface is assumed to be rigid, the boundary condition on 
the electron velocity is that the normal component is continuous (zero) at the interface. This is the 
condition used by earlier workers in the field [Yildiz and Silberg, 1964; Wait, 1964a, b; Johler, 1964; 
Field, 1956; Hessel, Marcuvitz, and Shmoys, 1962] . In addition, the usual conditions of continuous 
tangential electric and magnetic fields are used. 

In applying the boundary conditions to determine the unknown constants, it is immediately 
found that 

TT Tr ' 

which is to be expected from the discussion of equations (7). The equations governing the remain- 
ing terms can be written in the compact form 



an 



an 



a.u 



aii 



an 



«32 








~E Lr 

EoL 




bi 


023 




F' 

Eol 


= 


b 2 


aw 




F' 

Eol 




b. 



(32) 



fhere 



an = cos 6 + R sin d — iqfR sin #, 
an = sin <p\ + T cos <p\ — 7$, sin <pi , 
a2i = sin 9 — R cos 0, 
an =— cos <pi+T sin (pi, 

a23 = C08 <p 2 , 

a-.u =—t)lR, 
a32 = —r)Tu 

a33=—y]T2, 

6i =-R sin (9 + cos + t/£/? sin 0, 
bi =— sin 6 — R cos 0, 
bs =—r) L R. 



813 



In these equations, the indexes of refraction, t\l and tjt i, 2, have been introduced. 



7 w 
c 

j w 

Kt = — TJT' 
C 

The solution of (32) leads readily to the expression 

1f=T {2 cos 6 sin 6riTi(- l-R 2 ) + 2r)lr) T iR 2 sin 6 cos 6 

-\-2t)lR 2 cos 0[sin (fi + T cos </?i — t^ sin <pi] -\-2tjlR cos 0[— cos <pi + T sin <£i]}, (33a) 
where 
A= [cos d + R sin d — rjfR sin 0] [— T77^( — cos (£1 + T sin <£>i) + rjri cos ^2] 

+ [sin <p\ + T cos <£>i — 17^ sin <pi] [— 7]lR cos <p2 + ^^(sin 6 — R cos 0)] . (33b) 
In the case of zero static magnetic field, this reduces to 
E' TT _ — 2t)t\ cos 6 sin 



Eol 197-1 [cos 6 cos (f2~y]Ti sin sin <p\\ +sin sin <pi -hcos cos cpi 



(34) 



In (38), r)T2 has been set equal to unity. 

The conversion factors given by (33) and (34) represent the desired solution to the original 
problem. It is seen that conversion depends significantly on the angles. In the case of vanishing 
electron gyrofrequency, there is no radiation into the vacuum for normal incidence. This is not 
so for finite gyrofrequency, there being radiation into the plasma by virtue of the transverse wave 
which is generated in the plasma. 

An interesting aspect of the above development involves the direction of propagation. Every- 
thing so far developed holds equally well for complex angles. In fact, if the angle <p 2 is to be real, 
the angle of incidence, 0, must be complex. This fact comes out of Snell's law, 

Ul sin = kr\ sin <p\ — kri sin ip2. 

Since km is real and k L is complex (v ¥" 0), must be complex if <p 2 is to be real. This result is not 
surprising and in fact may be expected from earlier work. If is real, <p 2 is complex. The deci- 
sion as to how to choose the angles cannot be answered here. Instead consideration should be 
given to each of the two cases; i.e., real (<p2 complex) or complex (<p 2 real). Indeed, both cases 
can be expected to exist [Caron and Stewart, 1964] and, in fact, both angles may be complex. 
This is implicit in earlier published work [Wait, 1964b]. 

If consideration is given to the case where ip2 is complex, then primary concern should be di- 
rected toward the factor which determines the attenuation in the ^-direction. That factor will be 
dependent upon the incident angle as well as the propagation constants. The determining factor 
is the imaginary part of cos <p2 and the propagation constant in vacuum. By Snell's law, 



I £2 

k T 2 cos ip2 — ±k T 2 y 1 — jjr sm2 (35) 
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the choice of signs being made such that the field decays properly at large distances from the 
interface. In the present case, the lower sign is taken. 

It is interesting to note that in the absence of a static magnetic field and collisions, the term 
given by (35) is either real or imaginary, depending on the ratio of the plasma frequency to the 
radian frequency. In that case, 



exp (— ixkr2 cos ^2) = exp 



("*cVS 



1-- 



.-.). 



(36) 



If (Op < o>, then except for f-^J near unity and very small incidence angles, 0, the wave will decay 

exponentially with x because (c/u) is usually much greater than unity. 

There is an additional factor which lends insight to the nature of the conversion. This is 
the concept of power conversion. It is anticipated that this concept will be examined in a future 
publication. 



7. Discussion 

The variation of the electric field conversion factor, which describes the conversion of longi- 
tudinal (electroacoustic) field to electromagnetic field, is given in figures 2 through 6. The curves 
are parametric in electron density and are given as functions of frequency and (real) angle <p2. 
The value of - is taken to be 10~ 3 throughout (T is of the order of 10 4 deg Kelvin). The ordinate is 
a dimensionless conversion factor and can be related to perturbation pressure by a constant. 
Thus, if the incident longitudinal wave has a magnitude of 1 V/m, the associated pressure pertur- 
bation is 4.52 X 10" 7 |A:/jN/m 2 . The magnitude of k L varies from about 10" 3 at the lower fre- 
quencies to about 10~ 6 at the higher frequencies for electron densities of about 10 per cm 3 . At 
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FIGURE 2. Dimensionless conversion factor showing 
magnitude of electromagnetic field transmitted into 
vacuum for incident electroacoustic wave of unity 
magnitude. 

B o = 0; frequency =1 kc/s; collision frequency = 10 5 sec -1 . 
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Figure 3. Dimensionless conversion factor showing 
magnitude of electromagnetic field transmitted into 
vacuum for incident electroacoustic wave of unity 
magnitude. 

B o = 0; <^2 = 80°; collision frequency = 10 5 sec -1 . 
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Figure 4. Dimensionless conversion factor showing mag- 
nitude of electromagnetic field transmitted into vacuum 
for incident electroacoustic wave of unity magnitude. 

#o = 0.5 Gauss; frequency = 1 kc/s; collision frequency = 10 5 sec -1 . 
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Figure 5. Dimensionless conversion factor showing mag- 
nitude of electromagnetic field transmitted into vacuum 
for incident electroacoustic wave of unity magnitude. 

B — 0.5 Gauss; <^ 2 = 80°; collision frequency = 10 5 sec -1 . 
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FIGURE 6. Longitudinal-transverse coupling factor as a 
function of frequency for various electron densities; 
B = 0.5 Gauss; <p2 = 80°; collision frequency = 10 5 sec -1 . 

1 kc/s, |&l| has a value of about 10~ 4 . The results are given for the case of zero magnetostatic 
field (figs. 2-3) and also (figs. 4-5) for the case of a magnetostatic field having a magnitude of 0.5 
X 10" 4 Wb/m 2 (0.5 Gauss). The latter curves should be supplemented with the curve which shows 
the magnitude of /?, which is the magnitude of the transverse component of the quasi-longitudinal 
wave present in the plasma, for unity longitudinal component, by virtue of the static magnetic 
field, figure 6. Thus, some of the electromagnetic energy transmitted into the vacuum is generated 
by the transverse component incident on the boundary. 
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The dips which can be seen in the field conversion factors occur approximately at the plasma 
frequency in each case without the magnetostatic field. The introduction of the magnetic field 
causes this dip to shift slightly toward a lower frequency, especially for the lower electron densities. 

A problem complementary to that treated in section 3 concerns propagation in the direction 
of the static magnetic field. In that case, equations analogous to (7) can be found which show 
that there is no thermal correction to the index of refraction for the transverse components of the 
wave. It is then the longitudinal component which may propagate independently, the transverse 
components being coupled through the static magnetic field. 

8. Conclusions 

The E field conversion from an electroacoustic to an electromagnetic wave has been shown 
to be generally an increasing function of frequency and angle of transmission, up to one Mc/s. 
There is a slight dip at approximately the plasma frequency and this dip moves down the frequency 
scale when the magnetic field is introduced. The magnetic field tends to enhance the coupling 
at the small angles of transmission but detracts at the grazing angles. 
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